RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 69D, No. 2, February 1965 



On the Use of Refractive Index Diagrams for Source- 
Excited Anisotropic Regions* 

Leopold B. Felsen 

Contribution From the Department of Electrophysics, Polytechnic Institute of Brooklyn, Brooklyn, N.Y. 

(Received September 8, 1964) 

The use of refractive index diagrams in the study of plane wave propagation in un- 
bounded and layered anisotropic media is reviewed and then extended to account for certain 
aspects of the radiation from confined source distributions. Such items as the radiation 
condition, saddle point location, focusing effects, lateral ray trajectories, Cerenkov radiation, 
and others are interpreted via the refractive index plots. While some of this material is 
available in various technical publications, the aim here is at a more unified presentation. 



1. Introduction 

Since the problem of propagation of plane waves 
in a wholly or piece wise homogeneous anisotropic 
medium already involves considerable complexity, 
it may be anticipated that the difficulties are com- 
pounded further when spatially confined sources are 
present. The thorough understanding of the be- 
havior of plane waves is an essential prerequisite 
since the fields caused by localized distributions of 
sources may be synthesized by plane wave super- 
position. While the entire wave spectrum is required 
for a detailed description at any point in space, the 
field at observation points "far" from the source 
region, where a few propagating plane waves suffice 
to describe the essential features of the local field 
structure, may be characterized in a much simpler 
form. Jn analytical terms, the distant field is found 
from an asymptotic evaluation of the exact plane 
wave integral representation, the major contribution 
to which arises from the stationary points in the 
integrand. Each of the stationary points selects a 
particular plane wave which then appears in the 
asymptotic approximation. It may be surmised 
from these considerations that the selection of the 
appropriate waves appearing in the far field is a 
matter of basic concern. 

While these observations apply also to the iso- 
tropic case, their resolution is trivial in this instance. 
Since all propagating plane waves ~exp(ik-r) are 
characterized by the same wave number k = k n, 
where k is the wave number in vacuum and n is the 
(constant) refractive index in the medium, the field 
at a distant observation point P is furnished by that 
wave whose wave vector k is parallel to the radius 
vector from the source region to P. This situation 
no longer obtains in an anisotropic medium where 
k depends on the direction of propagation ; it is also 
necessary to distinguish between the directions of 
phase propagation (given by k) and power flow 
since these do not in general coincide. There arise 
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then some basic problems even about the deter- 
mination of the contributing plane waves in the far 
field, apart from the specification of their amplitude 
and polarization. 

Owing to the generally complicated functional de- 
pendence of k or n on the propagation angle [for the 
magneto-ionic medium, see Budden, 1961], it has 
been found useful to construct plots of these quan- 
tities which are termed wave vector and refractive 
index plots, respectively. From these diagrams, 
which have the form of multibranched surfaces, it 
is possible to ascertain the value of k and also the 
corresponding direction of power flow for any given 
spatial orientation. While refractive index surfaces 
have been employed for some time in connection 
with plane wave propagation, their use for radiation 
problems is of relatively recent origin. The latter 
aspect forms the motivation for this paper which 
begins with a review of some basic notions concern- 
ing the wave vector, group velocity and Poynting 
vector, and then proceeds to the study of radiation 
problems in an unbounded, homogeneous, lossless 
region. It is shown how the refractive index plot- 
may be utilized for the imposition of the radiation 
condition on the exact integral representation; the 
determination of the stationary points in the inte- 
grand, and therefore the selection of the plane wave 
constituents in the far field as well as the delineation 
of their domains of existence (some waves may appear 
only in limited sections of space) ; and the specifica- 
tion of those spatial regions wherein field enhance- 
ments (focusing) may take place (these effects are 
not observed in isotropic configurations) . The source 
configurations include point sources, line sources with 
and without progressive phase variation, highly di- 
rective distributions, and charges in uniform motion 
(Cerenkov effect). In all of these instances, the re- 
fractive index diagram furnishes in a simple and 
direct manner information about the frequently tor- 
tuous ray structure (trajectories of power flow) i\\h\ 
grants an insight into the radiation mechanism (but 
not directly into the polarization or angular intensity 
variation of a given ray species). These considera- 
tions are then extended to account for the presence 
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of a plane interface and for such phenomena as non- 
specular reflections, focusing, and lateral waves. To 
avoid the use of usually complicated analytical for- 
mulas in justification of results deduced for aniso- 
tropic configurations, a reasonably detailed applica- 
tion of pertinent concepts is given first for the 
thoroughly familiar isotropic case wherein the sur- 
faces have an especially simple shape. The extra- 
polation to the anisotropic regime should then be 
fairly evident. It should be mentioned that many 
of the aspects entering into the discussion have been 
noted variously in the literature to which reference 
is made; the principal aim in the present writing is 
at a unified presentation. Finally, while specific 
mention is made only of the magneto-ionic problem 
(single-species, zero temperature plasma), the same 
considerations apply also more generally [Light hill, 
I960]. 

2. Refractive Index Surface, Wave Vector, 
and Ray 

A space-time dependent electromagnetic field in 
an unbounded homogeneous medium may be repre- 
sented as a continuous superposition of plane waves 
which are characterized by a frequency w and by a 
wave vector k=x k x + j k 1/ -\-z k z . If the medium 1 
may be characterized at a given frequency in terms 
of a tensor permittivity e and a scalar permeability 

Mo (this applies, for example, to the magneto-ionic 
situation), each time-harmonic plane wave constitu- 
ent descriptive of the electric field, 

E(r) = Ae« k •'-**>, (1) 

must satisfy the vector wave equation, 



VXVXE-W-E=0, 
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(2) 



A is a constant vector, c the propagation speed in 
vacuum, e' the normalized permittivity, and r= 

x o^+yo2/+z 2 the position vector. The magnetic 
field may be derived in terms of E from the Maxwell 
equations. The parameters k and co in (1) are not 
independent but are connected by the dispersion 
relation obtained upon substituting (1) into (2) 
(note: v— nk): 



det Q=0, 



Q=kX(kXl)+^€' 



(3) 



where the unit dyadic 1 is defined so that 1 • A=A • 1 = 

A. Since (3) makes possible the solution for one of 
the parameters k x , k y , k z , co in terms of the other 
three, there exist two essentially different ways of 
representing this dependence: 

(a) w=w(k z , k yj k z ), or (b) k z =k z (k x , k y , w) . (4) 



i To simplify th^ discussion, the medium is assumed to be lossless. 



In case (a), the frequency variable is regarded as a 
function of the wave vector k whereas in (b), one 
of the wave numbers, k z , is a function of the remain- 
ing wave numbers k x , k y , and of to. The representa- 

A 

tion of the space-time dependent field E(r, t) via case 
(a) involves a triple Fourier integral in k-space 
[Ginzburg, 1960; Stix, 1962] 

E(M) = r dk x f" dk v P eft J(k)«* k •'-•«*, (5a) 

%) — 00 %J — oo c/ — CO 

while the one in case (b) leads to a representation in 
terms of "guided waves" along z [Budden, 1961], 
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with k t =k— Zo&g, p =r — z z . These integrations must 
be carried out separately for all possible solutions 
w(k) or k z (k t , w). 

The representations (5a) and (5b) are useful for 
different purposes. Energy transport properties of 
a quasi-monochromatic signal, expressed in terms 
of the group velocity vector \ g , are readily deduced 
from (5a). In this instance, A(k) has a sharp peak 
at a wave number k' corresponding to co =co(Z; ), and 

A 

the signal E(r, t) is a maximum along trajectories 
defined by the stationary phase points of the inte- 
grand for which the individual plane waves inter- 
fere constructively: 



V^(k)=0 at &', ?(k)=k • r-w(k)*, 



(6) 



with V Jt sXb(d/dfc)+y (d/dfc lf )+Zo(&/dJfc f ) denoting the 
gradient operator in &-space. With y g = r/t along 
the trajectory, one finds from (6) that 



▼,=V*w. 



(7) 



The direction of energy flow, parallel to \ g , is called 
the "ray direction" and differs in general in an 
anisotropic medium from that of the wave vector 
k which is perpendicular to the equiphase surface 
<p= constant. It is in fact evident from (7) that 
v^ is perpendicular to the surface co(k)= constant, 
i.e., to the constant frequency surface descriptive 
of the solution of the dispersion equation (3). Via 
(4b), the equation k 2 =k z (k x , k y ; a?) describes for a 
fixed co a surface which contains the endpoints of 
the wave vector k=\k\, and this plot in &-space is 
normalized conveniently by introduction of the 
refractive index n(k x , k y ; w): 



7 w 
k=—n. 
c 



(8) 



Apart from the normalization factor k =cx)/c, the 
wave vector and refractive index diagrams evidently 
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contain the same inform at ion at a specified frequency. 
Since the group velocity vector \ g (or the ray vector) 
for an essentially monochromatic plane wave has 
been shown to be perpendicular to the refractive 
index surface, it follows that k and v g are parallel 
everywhere only when |k| is independent of direc- 
tion, i.e., when the surface is a sphere. This obtains 
in an isotropic, but not in an anisotropic, region. In 
what follows, we shall be concerned almost exclu- 
sively with solutions of the dispersion equation 
which yield real values of k x , k y , k z corresponding to 
propagating plane waves in (1) ; while nonreal k 2 may 
arise in (5b) in view of the infinite range of k T , k v 
(see sec. 3.1), these evanescent fields do not carry 
energy away from the source region. It may be 
noted that monochromatic processes are described 
conveniently in terms of (5b) since the co-integration 
is then absent (the co-dependence of A(k t , co) is in 
the form 5(co— co ), where co is the operating fre- 
quency). Even these simple considerations illus- 
trate the utility of one or the other formulations in 
(4) and (5). 

For an explicit evaluation of the group velocity 
in terms of the refractive index n(k; co) , it is convenien t 
to proceed from the implicit differentiation formula 






i=x, y, z, 



(9) 



applied to f(k; u)=k—nu/c=0, which may then be 
solved for the components (dco/d/:,) of v f/ . If the 
^-surface is rotationally symmetric about one of the 
axes, say k z (this corresponds in the magneto-ionic 
case to an applied steady magnetic field parallel to 

z), it is convenient to define the polar angle d via 

sin ~d = k p /k, cos 6=kjk, and one may verify without 
difficulty that [Ginzburg, 1960; Budden, 1961; Stix, 

1962] 



cos a^— (noo) 



, v gk =r g cos a, v g d=v g sin a, 



where 



tan a= 



ndd 



(9a) 



(9b) 



In these expressions, //(#; co) is regarded as a function 
of d and co, v gk and v g e are the components of y g along 
and perpendicular to the wave vector k, respectively, 
and a is the angle between y g and k. One may show 
that |a|<C7r/2 — an observation important for the 
imposition of a radiation condition for problems 
involving confined source distributions in a magneto- 
ionic medium. Equation (9b) confirms that \ g is 
perpendicular to the n(d) surface. 

It is also of interest to state a simple relation 
between the group velocity v g and the time-averaged 
Poyn ting vector S = Re (EXH*) for a monochromatic 
plane wave propagating in a lossless, dispersive, 



anisotropic medium: 



S=Wv e , 



(10) 



where W is the average stored energy density, 



W- 



■H--S' 



coe(co)] • E+/z 



H| 2 } 



(10a) 



While certain aspects of (10) for the magneto-ionic 
medium may be verified by direct but tedious 
calculation [Hines, 1951 ; Abraham, 1953], its validity 
may be confirmed from a more basic and elegant 
analysis [Stix, 1962], the details of which are not 
presented here. The importance of (10) stems from 
the recognition that the ray direction in a lossless 
medium may be calculated either from S or \ g , with 
the former often the more convenient once the field 
evaluation has been completed. Evidently, the 

vector S is also perpendicular to the n(d) surface. 

It should be clear from the preceding discussion 
that the refractive index diagram contains pertinent 
information required for the determination of the 
energy flow characteristics in a monochromatic plane 
wave characterized by a given wave vector k. \\ hile 
the considerations so far have been restricted to an 
unbounded homogeneous medium, the refractive 
index plots may also be utilized effectively to 
ascertain the directions of the reflected and refracted 
rays at a plane interface between two media with 
different properties [Budden, 1961; see also sec. 5.3], 
By partitioning a slowly varying medium into a 
sequence of homogeneous layers, one may in this 
manner chart the progress of a continuously refracted 
ray — a procedure which lias found application in 
ionospheric propagation theory. It is the purpose 
in the remainder of this presentation to show how the 
refractive index plots may be used to advantage in 
predicting salient radiation characteristics of confined 
source distributions in unbounded and bounded 
anisotropic regions. As mentioned in the intro- 
duction, the pertinence of the preceding analysis 
stems from the fact that the fields due to spatially 
confined excitation, when observed at some distance, 
behave locally like plane waves. It may be noted 
that other plots (for example, of the ray velocity 
c(n cos a) -1 , or the ray refractive index (?i cos a)) 
may be employed to schematize wave propagation 
characteristics in an anisotropic medium. However, 
the refractive index plot is the most pertinent to the 
analysis of integrals as in (5b) since it utilizes the 
yt-space directly. 

3. Spatially Localized Source Distributions 
in Unbounded Anisotropic Regions 

The relevance of the foregoing remarks to radiat ion 
problems may be explored by referring to the integral 
(5b) which, upon omission of the co-integration, 
represents in terms of a superposition of the previously 
described plane waves the response E(r) due to non- 
moving, arbitrarily prescribed time-harmonic sources 
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(the time factor exp (—itd) has been suppressed): 

E(r)= r dk x f " dkMkx,h)e t[k ^ +k y v+k ^rVK 

(11) 

It is recalled that a separate integration is required 
for each of the relevant possible solutions k 2 (k x , k y ) 
of the dispersion equation. 




AXIS OF 
SYMMETRY 
(d-C MAGNETIC 
FIELD) 



PLANE OF SYMMETRY 



Figure 1. Typical wave vector plot in k y ,k z -plane 

Ti, 2— points of inflection. 

Dark portion— satisfies radiation condition for 2>0. 

Light portion— satisfies radiation condition for z<0. 



3.1. Radiation Condition 

The function k 2 (and A which generally depends 
on k z ) is multivalued in the k x , k y variables so that 
its definition is essential for the unique specification 
of the radiation integral (11). If the refractive 
index surface (or the wave vector surface k = k n) is 
of finite extent (closed surface), k z will be nonreal 
for sufficiently large (k x , k y ). Let us suppose that 
the excitation is in the form of a point source lo- 
cated at r=0, in terms of which any distributed 
driving function may be synthesized. Convergence 
of the integral then requires that Im k 2 ^0 for z^O. 
Physically, this radiation condition implies that the 
nonpropagating wave fields decay away from the 
source region. If the medium contains some losses 
so that k z has an imaginary part for all (k x ,k y ), the 
above-stated condition suffices to define k z every- 
where and assures the vanishing of the field at r-y»« . 
In analytical terms, the branch point singularities 
of k z are then displaced from the real k x ,k y axes, 
thereby making the integrand unambiguous along 
the entire integration path. 

To effect a definition of k z when this function is 
real, one may approach the limit of vanishing dissi- 
pation from the above-mentioned lossy case, and 
the required distortion of the integration path as 
the branch point singularities of k z approach the real 
axis is then made evident. This distortion of the 
path is equivalent to effecting the consistent defini- 
tion of the multivalued function k z for complex and 
real values. Alternatively, and frequently more 
conveniently, the proper definition of k zy when real, 



may be achieved directly from the refractive index 
plot [Arbel and Felsen, 1963]. Since the radiation 
condition requires that each constituent propagating 
plane wave transport energy away from the source 
region, 2 the permissible portion of the refractive 
index surface for z^>0 is that for which the surface 
normals have a component in the +z-direction; it 
is important to recall in this connection that the 
sense of the normal is such as to make its angle with 
the wave vector k less than 90° (see fig. 1 , where a 
section through the k y ,k z plane is shown). The 
converse argument applies when the observation 
point lies in 2<0. Having thus identified the 
appropriate branches of the refractive index plot, 
one may determine by inspection the correct alge- 
braic sign of k 2 . As noted from figure 1, it may 
happen that &/C0 when 2f>0 (rays S or S 3 ), thereby 
leading to a phase progression along { — z) while 
energy advances along (-\-z). This corresponds to 
a "backward wave" with respect to the 0-axis, a 
common occurrence in anisotropic media [see also 
Clemmow, 1963]. 

If the refractive index surface has an open branch 
so that k 2 is real for all values of k x and (or) k y , the 
preceding arguments may be applied as well. The 
integral (11) is then no longer exponentially con- 
vergent, and the field may grow without limit in 
certain spatial directions. This aspect is explored 
further in section 3.3b. 



3.2. Saddle Point Condition 

To effect a reduction of the integral (11), it is 
noted that the far fields are contributed essentially 
by those area elements in the k x , & y -plane surrounding 
the stationary points in the exponent. The station- 
ary phase condition 



x__5^ £ y_ dk 2 

z dk x z dk y 



(12) 



selects those values of k x and k y for which the normal 
to the refractive index surface points in the direction 
from the origin in the source region to the distant 
observation point, r. Since the normal to the 
surface has been seen to specify the direction of 
power flow, one recognizes that the major contri- 
bution to the far field arises precisely from those 
plane waves (rays) which carry energy from the 
source to the observation point along a straight line 
trajectory. Conversely, the refractive index plot 
may be utilized to determine the real stationary 
(saddle) points in the k x ,k y wave number space: 
one finds all points on the surface having a normal 
parallel to the prescribed radius vector r (see for 
example the four rays S 2 . . . S 4 in fig. 1), and 
reads off the corresponding values of k x ,k y which. 



2 The radiation condition, which requires the outward flow of total energy, 
must be satisfied for arbitrary source configurations, i.e., arbitrary A in (11). 
One then concludes that each plane wave constituent must individually satisfy 
this condition [see Arbel and Felsen, 1963]. 
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locate the saddle points, remembering the previously 
mentioned radiation condition that the angle be- 
tween k and the ray is less than 90°. These graphical 
considerations therefore aid substantially in the 
visualization of the saddle point configuration 
pertaining to the integral (11) [Lighthill, 1960; 
Arbel and Felsen, 1963; Mittra and Deschamps, 
1963; Kogelnik and Motz, 1963]. 

3.3. Ray Amplitudes 
a. Curvature of Refractive Index Surface 

The refractive index plots furnish information 
not only about the multiplicity and propagation 
characteristics of the rays reaching a distant observa- 
tion point, but also about the ray amplitudes. 
The asymptotic evaluation of the double integral 
(11) may be carried out by adding the contributions 
from the vicinity of all appropriate stationary points 
satisfying condition (12). Since the choice of 
coordinate axes with respect to the dispersion surface 
has been left arbitrary, we may select an orientation 
so that the 2-axis is parallel to the radius vector r; 
i.e., x=y=() in (11), and the saddle points are 
defined by those points on the surface with normal 
parallel to z. In the vicinity of a saddle point 
(k xj ,k yj )=k tj , the presumedly slowly varying func- 
tion A may be approximated by A(k xj ,k vj )=Aj, and 
the function k z in the exponential may be repre- 
sented by the first few terms in its power series 
expansion about (k xj) k VJ ). In view of (12), with 
y=z = 0, the linear terms in (k x —k xj )=% and 
(k u — k yj )=r) are absent. By orienting the {k x ,k v ) 
coordinates so that they coincide with the principal 
directions of curvature of the surface at the saddle 
point, one eliminates the i^-term and may write to 
terms of second order, 

k z (k x , k y )=k 2j + hK xj ? + hK y tf+ . . ., (13) 

with K xj =d 2 k z /dk 2 x ] j and Kyj=d 2 k z /dkl]j representing 
the associated curvatures, and k zj =k z (k xj , k yj ) . Thus, 
the asymptotic approximation of E(r) for large 
values of z is given by [Lighthill, I960]: 

E(r)~S Xje ik ^ z f e iaK *i?*dZ f e izK v^ l2 dy] (14a) 

;' J — oo J — oo 



Jh y- A , «,,. e'f **+*■**■'» (l_ 



(I) (14b) 



where Kj = K xj K yj is the Gaussian curvature at the 
jth contributing saddle point included in the sum. 
The detailed structure of Aj depends both on the 
medium parameters (i.e., the refractive index profile) 
and on the nature of the source configuration. 
While the preceding considerations yield all the 
possible rays which may propagate from the source 
to the observation point, it should be kept in mind 
that specific source configurations may excite only 
some of these. 

Equation (14b) states that the distant field 
radiated by a confined source distribution is com- 



prised of radial rays along which the field decays 
according to (1/r). The angular pattern (i.e., the 
amplitude) of the field along a given ray is governed 
by the coefficient A ; which is generally well behaved, 
and by the factor Kj 1 ' 2 . Disregarding A ; , the 
amplitude evidently increases when the ray originates 
from a "slowly curved" segment of the dispersion 
surface, and ic appears to grow without bound 
when one of the principal radii of curvature vanishes. 
A closer investigation of the refractive index plot 
reveals that Kj-^0 implies the coalescence of two 
saddle points and therefore the strong interaction 
of two almost identical ray species (see fig. 1, where 
ky\-^k y2 when the radius vector r from the source to 
the observation point becomes perpendicular to the 
curve at the inflection point T Xj similar considera- 
tions apply to T 2 ). This situation requires a modi- 
fication of the asymptotic evaluation of (11) wherein 
it was assumed that the saddle points are simple 
(K xj ,K vj 5*0) and separated. The calculation when 
Kj—>0 leads to a radial dependence involving Airy 
functions and to field strengths enhanced over the 
(1/r) dependence; details have been given elsewhere 
[Arbel and Felsen, 1963] and will not be repeated 
here. We should merely like to emphasize in the 
context of the present discussion that the occurrence 
of field enhancement by virtue of the interaction of 
rays arising from two adjacent stationary points 
may be predicted directly from the presence of 
points of inflection in the refractive index diagram. 

b. Open Branches of Refractive Index Surface 

If the refractive index surface contains unbounded 
branches, there will exist certain ray directions 
corresponding to saddle points k xj ^>co and (or) 
kvj~ >0 ° • Since k z is now real over the entire integra- 
tion interval, the integral (11) is no longer exponen- 
tially convergent, and E(r) may actually diverge 
when the saddle point moves to infinity. The 
nature of the singularity depends on the source 
configuration; it is most pronounced for point 
source excitation but may actually be made to dis- 
appear for sufficiently smooth source distributions. 
Details of the analysis are available elsewhere 
[Arbel and Felsen, 1963; Felsen, 1964a] and it 
suffices for the present to note that the singular ray 
directions may be inferred easily from the refractive 
index diagram (see ray S' in fig. 1, as &— ><» along the 
open branch ; the singular ray direction is perpendic- 
ular to the asymptote). 

3.4. Domain of Existence of a Particular Ray Species 

If the refractive index surface possesses undulations 
("bumps" or "dimples") or unbounded branches, 
not all regions of space are illuminated by the same 
ray species. For example, if X and 6 2 denote the 
angles between the positive & 2 -axis and the normals 
at T x and T 2 , respectively, in figure 1, ray Si reaches 
those observation points in the first quadrant of 
the y, z plane whose angular deviation 6 from the 
2-axis lies in the interval 0<6<6i. Ray S 2 is con- 
fined to 6 2 <d<6i, ray S 3 exists in 6 2 <6<t/2, and ray 
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S 4 in dz<d<w/2 with 3 representing the inclination 
of the normal to the asymptote of the open branch. 
The angular domain wherein a ray species exists 
is its zone of illumination and the remaining region 
is the shadow zone wherein the fields on the ray in 
question are exponentially small. It is evident that 
the zones of illumination and shadow may be inferred 
directly from the refractive index plot from which 
it is seen, for example, that rays S x and S 2 are absent 
when 6i <6<tt/2 (the corresponding saddle points are 
complex in this range) while the shadow region for 
rays S 2 and S 3 is O<0<# 2 (see also fig. 8). 

4. Source Problems in Isotropic Regions 

To illustrate specific application of the preceding 
concepts in their simplest form, a number of radiation 
problems in isotropic configurations are reviewed 
and interpreted via the refractive index plots. In 
a homogeneous isotropic medium, ?=le, and one 
obtains from (3) the dispersion relation 



k=-^=k n, 
c 



k=^kl+Ic 2 y +Jc i z =^k 2 +k 2 , 



where n may depend on co but not on 6. As noted 
previously, the refractive index plot reduces to a 
sphere in this instance so that the directions of the 
wave vector k and the ray S coincide. 

4.1. Sources in an Unbounded Region 

a. Point Source 

The fields radiated by an arbitrarily oriented 
electric or magnetic current element in an infinite, 
homogeneous, isotropic dielectric may be derived 
by suitable vector operations from the scalar Green's 
function 

n ikr 



which lias the integral representation 



(15a) 



G-- 



I 8^J_„ d **J_ 



dk„ 



h 



k z =T>kl71?-kl-kl. (15b) 




Figure 2. Wave vector surface for isotropic medium: k = k n. 
Saddle point— k y \=k sin 0. 



(The above-mentioned vector operations account for 
the detailed structure of A in (11) for different 
source configurations.) While the solution for this 
trivial example is known in the closed form (15a), 
the integral (15b) shall be investigated per se for 
purposes of illustration. 

Because of the spherical symmetry, the coordinate 
system may be chosen so that x=0, thereby leading 
to the diagram in figure 2. In view of the dependence 
on \z\ in (15b), only the region 2f>0 need be con- 
sidered. k z must be positive imaginary when non- 
real, and it is evident from the diagram in figure 2 
that to satisfy the radiation condition, k z is positive 
when real. It follows from this definition that the 
integration path is indented into the lower and upper 
halves of the k x plane to avoid the branch points at 

k x =-yjkln 2 —k 2 y and at k x =—^k 2 n 2 —k 2 y , respectively, 
when k 2 y <Ckln 2 . Analogous considerations apply 
to the ^-integration. The same conclusions are 
reached when dissipation is assumed initially (Im 
t?T>0) and the lossless case is approached subse- 
quently. One observes, incidentally, that the branch 
points correspond to those portions of the refractive 
index surface for which the associated rays are 
confined to the z=0 plane (separation between rays 
going into z^>0 and 2<C0, respectively); this latter 
observation holds also under more general circum- 
stances where these singularities are no longer 
defined by the simple condition k z =0. 

From the discussion in section 3.2 and figure 
2, it follows that the only saddle point of the inte- 
grand in (15b) (with x=0) is located at k yl = k 
sin 6, where Q=tsm~ 1 (y/z) specifies the angular 
location of the observation point. While this 
conclusion is reached at once from the geometrical 
construction, it may of course be verified from (12) 
by direct calculation. Each distant observation 
] joint is therefore reached by a single radial ray 
with wave vector k. 

Since the principal curvatures of the spherical 
surface are constant and equal to (— 1/&), one 
obtains directly from (14b) that the asymptotic 
approximation to G is given by the expression 
(15a) (which in this instance happens to hold for 
all r). None of the complicating features due to 
points of inflection or open branches of the re- 
fractive index surface arise here. It should be 
recalled that the angular variation of the ray ampli- 
tude also depends on the nature of the source 
and emerges from the derivation of the field via 
the scalar Green's function G. 

b. Line Source 

Constant phase. If the excitation is in the form 
of electric and (or) magnetic currents distributed 
uniformly along a line, say the 7/-axis, then the 
entire field structure is independent of the ^/-coordi- 
nate, and the constituent plane waves descriptive 
of this field configuration must be characterized 
by k y =0. Consequently, the effective portion of 
the refractive index diagram is the curve formed 
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by the intersection of the sphere in figure 2 with 
i he plane ^ = 0. The corresponding surface normals 
are confined to the x, z -plane so that energy leaves 
the source along radial rays which lie in a plane 
perpendicular to the source axis. The polarization 
of the field, or its angular intensity distribution, 
depends, of course, on the details of the source 
configuration. 

Linearly progressing phase. If there is impressed 
on the previously mentioned line distribution a 
phase variation exp (ik y] y), where k, n is a real con- 
stant, then the resulting fields in this ^/-invariant 
structure must exhibit the same ^-dependence. The 
constituent plane waves are therefore characterized 
by ky = k yh which requirement reduces the double 
integral in (11) to a single one over k x (the in- 
dependence of A(k x , k y ) is in the form b{ky—k lA )). 
1 lie radiation characteristics are now inferred from 
the intersection of the plane k y =k yl with the refrac- 
tive index surface in figure 2. One observes that 
if l&yil^i, energy emanates from the source along 
rays which lie on a circular cone Inning an apex 
angle ^=cos _1 (k yl /k) with respect to the source axis 
(fig. 3). It is evident from figure 2 that real rays 
exist only when \k yl \<ik. For \k yl \y>k, there is no 
real solution for k z and radiation does not occur. 
The intervals \k n \<ik and \k y i\^>k distinguish phase 
velocities along the source distribution which are 
larger and smaller, respectively, than the velocity of 
light in the medium; as is well known, the former 
wave types radiate whereas the energy in the latter 
is bound to the source region. 

The validity of the preceding remarks is easily 
verified from the scalar Green's function 



Q* e * vl V H U y(k 2 -k 2 yl ) (X 2 +Z 2 )} 



(16a) 



=^ e**i* J -j— dk XJ h=jk 2 -k 2 x -k 2 yl , 

(16b) 

in terms of which the electromagnetic fields may be 
derived. For example, if the source is a line of elec- 
tric currents of unit strength, then the non vanishing 
field components are 



E ' = h(to? + b?) G > ER= ^Mdy' H * 



(16c) 



with R=^x 2 +z 2 and </>=tan -1 (z/x) denoting polar 
coordinates in the plane transverse to y. Using the 
asymptotic approximation of 0, one finds E^ (vX<l>o)k 
(k 2 —Ey X )~ ll2 E, h H^-\/e/iJLvxE, where v is a unit 
vector along one of the rays sketched in figure 3. 
The Poynting vector is therefore parallel to v. For 
k !n = 0, one recovers the case described in the first 
paragraph. 




Figure 3. Ray configuration for phased line soura . 



c. Twc-Dimensional Source Distribution 

If the sources are distributed continuously over 
the entire .r,7/-plane with impressed phase variation 
exp (ikxix) exp (ik yl y) , where k xl and k yi are real con- 
stants, then the field is described by a single plane 
wave with (k xU k yl ). The direction of energy flow 
is again inferred from that of the normal to the 
refractive index surface at this point. It is noted 
that radiation takes place only if k 2 xl -\-k 2 y i<^k 2 . 

d. Highly Directive, but Confined, Source Distributions 

The preceding phased one- and two-dimensional 
source distributions of infinite extent are idealiza- 
tions which may, however, serve as approximations to 
large phased arrays or apertures of finite extent. 
Consider, for example, a progressively phased line 
source as in figure 3, confined to the interval —a<C 
y<Ca. The amplitude function A(k x , k y ) in (1 1) now 
contains the factor [sin {k y — k n )a]{k y — k y y)~ l } which 
has a sharp peak at k y =k y] if the source region is 
many wavelengths long. While this behavior is less 
violent than the 8(k y —k y i) corresponding to a->oo ? it 
nevertheless selects the plane wares with k y =k yl as 
those contributing to the major peak(s) in the 
radiation field. The previously determined rays 
therefore point in the direction of the major lobes in 
the radiation pattern so that the refractive index 
diagram may be employed to find the locations of the 
radiated beam maxima. 

e. Radiation From a Uniformly Moving Charge (Cerenkov Effect) 

If a point charge a moves with constant velocity v 
along the y-direction, the associated electric current 

A 

density J may be characterized as follows: 



i=J,qvb{x-x f )b{z-z f )b{y-vt). (17) 

The frequency spectrum J of this source distribution 
is obtained by taking the temporal Fourier transform 



s !f* 



e +iu, }dt=jo^S(x-x')d(z-z')e Uk o^", 



(18) 
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where k =o>/c and (}=v/c. This excitation for 
specified frequency corresponds precisely to the 
linearly phased line current discussed in section 4.1b, 
with amplitude (q/2ir) and impressed phase progres- 
sion k yl = k Q /(3. The radiation characteristics at any 
given frequency may therefore be inferred from the 
results in section 4.1b, and the space-time dependent 
field is then recovered upon carrying out the co-inte- 
gration in (5b) . 

Since the particle speed is always less than the 
speed of light in vacuum, i.e., /3< 1, one observes that 
k y i^>k . It is then obvious from figure 2 that no 
radiation will take place in vacuum for which the 
refractive index n=l. In fact, radiation is possible 
only when the charge moves in a medium with 
refractive index rf>(l//3) for which one may satisfy 
the condition kyi^k. If dispersion in the medium is 
neglected so that n is independent of frequency, then 
radiation at all frequencies will emerge along the 
well-known "Cerenkov angle" 

^cos-fl^cos-^l) (19) 

which may be determined graphically from the 
refractive index diagram (figs. 2 and 3). For a 
frequency dependent n, the Cerenkov angle will 
likewise be a function of frequency. 



4.2. Influence of a Plane Boundary 

In the preceding sections, the refractive index 
surface has been utilized for the determination of 
the power flow properties associated with various 
source distributions in unbounded media. The discus- 
sion is easily extended to accommodate the presence 
of a plane boundary which either terminates the 
region (for example, a perfect conductor) or separates 
it from another with different physical properties. 
The underlying analytical considerations are similar 
to those dealt with earlier. If it is assumed, for 
example, that the plane 2=0 separates two semi- 
infinite dielectrics with n=ni(z<^0) and n=n 2 {z^>Q) , 
then the fields of a ^-directed point source located at 
(0, 0, z'), z r< C0, may again be derived from a scalar 
Green's function G which has a different representa- 
tion in the two regions. In region 1 (s<0), 
contains a direct and a reflected portion [Stratton, 
1941] 



G=G[+G[', 



• fco /»oo „i(k T x+k v y±k x z) 

G[=^_ dk x j_Jk y e — 



G['= f dk x (" dk y T(k x ,k v ) 



e i(k z x+k y y-k 2l z) 

k zl 



. e ~^zl z 



(20) 
(20a) 
(20b) 



whereas in region 2 (s>0), the transmitted field is 



derived from 



G=G 2 = P dk x p dk y T(k x , k y )e' 



{kjjc+k y y+k z2 z) „ - ik zl zf 



k 2i ^kln 2 i-¥ x -kl. (21) 

Here, r and T represent the plane wave reflection 
and transmission coefficients, respectively, whose 
detailed form (for either E or H polarization) is 
readily evaluated but irrelevant to the present 
discussion. The important thing to be recognized 
is that the constituent plane waves in both regions 
are characterized by the same wave numbers k x , k v 
in the direction parallel to the interface — a condition 
required to assure continuity of the tangential field 
components across the boundary. The signs associ- 
ated with k zi and k z2 have been chosen in accord with 
the previously discussed radiation condition (sec. 
4.1); in (20a), they assure the outward power flow 
from the source plane z=z f , in (20b) the flow of re- 
flected power toward z= — &> y and in (21) the flow of 
transmitted power toward z= + ^. 

a. Incident, Reflected and Refracted Rays 

The familiar relation between the angles of inci- 
dence, reflection, and refraction of a plane wave, 
known as Snell's law, is inferred at once from the 
refractive index plots, a section of which is shown in 
figure 4. The propagating wave solutions corre- 
sponding to a given value k y = k yl (with ^=0 through 
proper choice of coordinates) may be ascertained 
from figure 4 wherein it is assumed that n{^>7i2. 
Evidently, four solutions are possible: two rays 
each progressing into 2>0 and 2<0, respectively. 
If a ray is incident on the interface from region 1 
(S[), the corresponding reflected ray is Si'; since 
the surface k=k x is rotationally symmetric about 
the & 2 -axis, the angles of incidence and reflection are 
equal (d{). Of the two ray solutions S 2 and S2' 
corresponding to the prescribed k yi in region 2, 
S^ must be discarded since it carries energy along 
the negative ^-direction and therefore violates the 
radiation condition associated with an incident ray 
from below. The relation between the angles of 
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Figure 4. Wave vector surfaces and rays for two different 
isotropic regions: ki = k ni. 
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refraction (6>) and incidence is obtained from the 
graph via the condition: k vX = k l sin 8 { = k 2 sin 2 , or 
sin &2=(ni/n 2 ) sin 6 X (SnelPs law). Since the two 
surfaces are spherical, the incident, reflected and re- 
fracted rays all lie in a common plane. One notes 
that the ^-branch is intersected only when k yl < Ck 2 ' ) 
for k yl ^>k 2j i.e., sin 0i>(n 2 /^i), k z2 is imaginary and 
no propagation takes place in region 2. The incident 
wave is then totally reflected. 

b. Lateral Rays 

When k yl =k 2 j the ray is incident at the critical 
angle 0; = sm -1 (n 2 /rii) and the refracted ray proceeds 
parallel to the interface (see right-hand portion of 
fig. 4). This parallel ray, called a "lateral" ray, 
may refract back into region 1, thereby providing a 
mechanism for energy transfer which is entirely 
different from that associated with a reflected ray 
[Brekhovskikh, I960]. In contrast to the discussion 
in section 4.2a which is relevant for an ordinary 
plane wave, the appearance of a lateral ray is con- 
nected intimately to the presence of a source of finite 
dimensions and owing to the continuous leakage of 
energy (see fig. 4), the amplitude of the field on a 
lateral ray is smaller than that on a reflected ray. 
While some further comments on tins ray species 
are made later on, it is to be emphasized at this point 
that the existence and trajectory of a lateral ray 
may be inferred directly from the refractive index 
plot. 

If the source is placed into the "thinner" dielectric 
(i.e., if n 2 >^i), the lateral ray is not excited since it 
is not possible to have 2 =7r/2 for any incidence 
angle in the range O<0i<7r/2. The intimately 
related phenomenon of total reflection is also absent 
in this instance. 

c. Asymptotic Evaluation of Reflected and Transmitted Fields 

While the direct field G[ may be evaluated as 
before (see (15a, b)), the reflected and transmitted 
fields require further attention. If we again set x=0 
for convenience, the saddle point in the k x integral is 
located at k xs =0. The saddle point condition applied 
to the remaining integration in (20b) yields: 



y+ 



bk 2l , | . bhi 



bk y 



*\+ 



dh 



1 3' 1=0 at Jc y . 



(22) 



Since (i>k zi /i>k y ) = — tan 6 U where 6 t is the angle 
of the normal to the refractive index diagram at 
{k Z i- k v ), the straight line through (y, z) defined by 
(22) may be interpreted geometrically as follows: 
The observation point (y f z) in region 1 is reached 
via the ray trajectories Si and Si' in such a 
manner that the angle of reflection is equal to the 
angle of incidence; this condition determines X and 
thereby the saddle point k ys via figure 4. The 
equality of the angles of incidence and reflection 
follows from the symmetry of the refractive 
index plot; the ray interpretation of the saddle 
point condition is unaltered for more general situ- 
ations where these angles may be different from one 



another (see sec. 5.3). Confirmation of the ray 
interpretation in figure 5 is also had from the expo- 
nential in (20b) evaluated at the saddle point: 
exp [i(ky S y+k zl \z\+k zl \z'\\. Since y=y { +y 2 , the phase 
change is precisely the one required by a plane wave 
to travel from the source to the observation point 
via the ray trajectory Si, Si'. 

Analogous considerations apply to the evaluation 
of the transmitted field via (21). The saddle point 
condition now reads 



y+ 



dk v 



dk z2 
~i)k v 



z=0 at k v , 



(2;?) 



and has a simple graphical interpretation: The 
observation point {y, z) in region 2 is reached via 
the ray trajectories Si and S 2 in such a manner that 
the relation between the angles 6 X and 6 2 is that spec- 
ified by the refractive index plot in figure 4. 

In summary, the saddle point determination re- 
quired for the asymptotic evaluation of the reflected 
and transmitted fields, and the subsequent ray- 
optical interpretation of the result, may be carried 
out with the aid of the refractive index diagram: 
The source and observation points are connected 
by rays which satisfy the plane wave reflection or 
refraction laws, and the corresponding values of 
k v yield the saddle points. From a knowledge of 
the saddle points, one may determine the ray ampli- 
tude via the remaining amplitude factor in the 
radiation integral. 

The saddle points associated with propagating 
wave solutions in both regions are restricted to the 
interval \k ys \ < Ck l and \k ll!i \<^k 2 for the reflected and 
transmitted fields, respectively. If k{^>k 2 , then the 
saddle points descriptive of the reflected field may 
cross the branch points at ±k 2 . Under these con- 
ditions, the asymptotic field contains another con- 
stituent arising from a branch point integration. 
This field type is the previously mentioned lateral 
wave whose existence and trajectory is readily pre- 
dictable from the refractive index diagram. As 
noted from figure 4, the lateral wave contributes 
only when d 1 ^>6[, thereby implying that the branch 
point contribution is absent when \/c ys \ < Ck 2 . This 
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Figure 5. Ray interpretation of saddle point condition for 
reflected and refracted fields (yi+y2=y). 
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a) original path 



ys 



b) steepest descent path 



Figure 6. Integration contours in complex k y -plane. 

a. Original path. 

1). Steepest, descent path. 



conclusion is confirmed from a study of the steepest 
descent path through the saddle point which involves 
a deformation of the integration contour away from 
the real axis, thus making unambiguous the domain 
of contribution from the branch point (see fig. 6). 

d. Other Source Distributions 

The preceding discussion has shown that the 
plane wave refractive index diagrams have direct 
application to the determination of the reflected or 
transmitted fields when a plane boundary is present. 
In essence, one first finds the ray structure in the 
unbounded region containing the source and then 
accounts for the interface by a ray-tracing procedure, 
with the trajectories of the reflected and transmitted 
rays ascertained from the diagram. With this in 
mind, the study of other source configurations is 
evident. For example, the highly directive sources 
described in section 4. Id excite strongly certain 
selected ray configurations whose characteristics in 
the presence of the interface are predicted from the 
refractive index plot. In this manner, one may 
find the directions of the reflected and transmitted 
beams. Cerenkov radiation due to a charge moving 
parallel to the boundary is also easily understood. 
If the charge moves in region 1 and ky=(k v //3) < Ck 2 
<C/:i, then radiation will occur in both regions at 
readily determined angles. For k{^> (k /P) >& 2 , prop- 
agation is possible only in region 1 and the energy 
is totally reflected at the interface, whereas for 
(k Q f/3)y>ki, no radiation takes place. If the charge 
moves in region 2, radiation escapes into both regions 
when (k /(3)<Ck 2 < Ck 1 , but the fields in region 2 are 
evanescent when &i>(& c /£0>&2. Since propagation 
in region 1 is possible, however, under the latter 
condition, the incident evanescent fields may couple 
to propagating waves at the boundary. 

4.3. Multiple and Gently Curved Interfaces 

If the region is comprised of homogeneous layers 
(an inhomogeneous region may be so approximated) , 
each layer has its own refractive index diagram and 
the previously described field matching procedure 
may be employed to chart the progress of a ray 



through the medium. The diagrams also reveal 
when a ray is "trapped" (ducting effect, bound 
waves): if a layer is bounded on both sides by 
media with smaller refractive index, it is possible to 
find wave solutions which propagate in the layer 
but decay outside. 

In view of the ray interpretation of the fields 
radiated by confined source distributions, the ray 
tracing method may be applied locally even when an 
interface is non planar provided that the curvature 
is small over a length interval equal to the local 
wavelength [Keller, 1958]. Under these circum- 
stances, the boundary at the point of impact of the 
ray is assumed plane, and the refractive index plots 
may then be used to determine the initial trajectories 
of the reflected and transmitted rays. 

5. Further Discussion for Anisotropic 
Regions 

The discussion in section 4 involving isotropic 
regions has served to illustrate the use of refractive 
index diagrams for the study of various familiar 
radiation problems. While the associated phenom- 
ena and their interpretation in terms of rays are in 
this simple case also easily deduced by other means, 
the refractive index plot offers distinct advantages 
in anisotropic regions where its contour may depart 
drastically from the spherical shape. With an 
understanding of the isotropic situation, one may 
proceed to apply the same considerations to the aniso- 
tropic case; the only difference is the substantially 
more complicated configurations of surfaces which 
depend on the form of the dielectric tensor *J in (3). 
For the magneto -ionic medium, they may be grouped 
into some eight categories [Clemmow and Mullaly, 
1955; Allis, 1961; Deschamps and Weeks, 1962] 
which range from two concentric ovals to the undu- 
lating and open-branched configuration shown in 
figure 1, as the applied, cyclotron and plasma fre- 
quencies (for electrons) take on various values. The 
variety of surface contours is further increased when 
additional species (for example, ions) are taken into 
account in the description of a magneto-plasma 
[Allis, Buchsbaum, and Bers, 1962]. For the present 
discussion, the specific nature of the surfaces is of no 
concern and it suffices to deal with a typical case (as, 
for example, in fig. 1). It should of course be evident 
that the multibranched character of the refractive 
index plots implies the existence of more rays than 
those observed in the isotropic case. Even in the 
magneto-ionic medium wherein only electrons are 
considered mobile, as many as four different propa- 
gating solutions for k z (k x , k y ) may be encountered. 
In view of the very complicated analytical expres- 
sion for k z in terms of k x and k y [Budden, 1961], it 
may be appreciated that much insight is gained from 
a plot of these quantities provided by the refractive 
index diagram. 

The problem of radiation from a point source in 
an unbounded anisotropic region has already been 
dealt with in section 3. The anisotropy exerts in- 
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teresting effects also on other source distributions 
which are readily explored. 

5.1. Line Source 

The power flow characteristics for a uniform line 
source directed along the y-nxis are inferred from the 
k v =0 section through the refractive index plot. If 
the y-axis coincides with the direction of the axis of 
symmetry (gyrotropic axis), the normals to the sur- 
face in the Jc y = plane are contained in this plane; 
the radiation therefore leaves the source in the radial 
direction and the rays are perpendicular to the line 
axis. For arbitrary inclination of the gyrotropic 
axis, however, the surface normals need not lie in 
the k y =0 plane. With reference to figure 1, for 
example, the two rays with k x =0, k y = point into 
the first (top ray) and third (bottom ray) quadrants, 
respectively, with intermediate positions occupied 
by rays with kx^O. Energy thus travels along 
trajectories which arc no longer perpendicular to the 
source axis although no phase variation is impressed 
along the source, a behavior in marked contrast to 
that observed earlier in the isotropic case. If (he 
k y = plane intercepts more than one branch of the 
refractive index surface, several ray species (with 
different orientation) may arise. 

These observations apply also when the source 
possesses a progressive phase variation ex p (//•„,//), in 
which instance the intersection with the plane 
k y =kyi is pertinent. For example, if k yX is so large 
that only the open branch of the surface in figure 1 
is intercepted, one obtains the ray structure shown 
in figure 7. In view of the requirement k-S>(), it 
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b) line source and rays 

Figure 7. Progressively phased line source inclined to 
gyrotropic axis {gyrotropic axis lies in y-z plane) . 

a. Wave vector surface and k y =k v i section. 

b. Line source and rays. 
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a) wave normal plot 



b) ray configuration 



Figure 8. Unphased line source perpendicular to 
gyrotropic axis. 

a. Wave vector plot. 

I). Ray configuration. 



is noted that the rays S 2 and S 2 in the k x =0 plane 
are "backward" in the sense that the direction of 
phase progression is toward the source axis whereas 
energy necessarily flows outward. Evidently, the 
rays are not confined to a circular cone as in the 
isotropic problem (fig. 3). 

Since 1 the refractive index diagrams for the mag- 
neto-ionic medium are rotationally symmetric about 
the gyrotropic axis, substantial simplification in the 
ray picture results when the line source is oriented 
along, or perpendicular to, the axis of symmetry. 
In the former case, already mentioned earlier, the 
ray configuration is also rotationally symmetric 
and the rays in figure 7, for example, lie on a cir- 
cular cone. In the latter instance, with k y =0, the 
effective portion of the diagram comprises the plane 
section in the k x , k 2 plane as illustrated in figure 8 
for the case of an open-branched surface. The rays 
are now confined to the x, z plane (i.e., they leave 
the source at right angles) and they illuminate only 
the shaded region in figure 8b). 

5.2. Other Source Distributions 

The preceding considerations concerning line dis- 
tributions (and the analogous surface distributions 
described in section 4.1c) are again applicable to 
spatially confined but highly directive sources (see 
section 4. Id). For example, a planar phased array 
in the x-y plane with k x =0, k y =k yh produces two 
narrow beams which point asymmetrically along 
Si and S 2 in figure 7b, whereas a linear phased 
array with k y =k y i produces beam maxima over 
a noncircular cone whose generators follow from 
figure 7 [Felsen and Rulf, 1963]. The result for 
the linearly phased line source (with k y ] = k //3 as in 
section 4. Id) may also be employed for the study 
of the radiation produced by a given spectral com- 
ponent at frequency co arising from a point charge 
moving uniformly along the y-axis. Because of the 
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highly dispersive character of the magneto -ionic 
medium, the radiation characteristics differ greatly 
for various frequency ranges. Once the refractive 
index plots are available, it is a simple matter to 
ascertain when radiation may take place (i.e., when 
the Cerenkov "coherence condition" may be satis- 
fied): one looks for possible intersections of the 
plane k y =(k o /0) with the wave normal surface and 
infers the corresponding ray angles [Clemmow, 1963; 
McKenzie, 1963]. 

It must be emphasized that these graphical 
methods for the determination of major features in 
the power radiation pattern apply directly only 
when a single ray reaches a given observation point. 
If several rays are involved, their combined power 
pattern (which may deviate from the individual 
behavior) must be determined. It may happen, 
however, that the given source configuration excites 
one ray species more strongly than the others, in 
which instance the pattern is determined essentially 
by the most strongly excited ray. 

5.3. Presence of a Plane Boundary 

a. Ray Interpretation of Fields 

If the gyrotropic medium is bounded by a plane 
interface at 2 = 0, the radiation integral (11) must 
be augmented by additional contributions which 
account for the reflected and transmitted parts 
as in (20) and (21). Since there will generally 
be several solutions for k 2 (k x ,k y ) as evidenced by 
the multi-branched refractive index plots, it is 
understood that a separate integration is required 
for each pertinent k z (corresponding, for example, 
to the ordinary and extraordinary modes). The 
resulting integrals are quite complicated [Barsukov, 
1959; Arbel and Felsen, 1963; Tyras, Ishimaru, 
and Swarm, 1963] but as in the isotropic case 
(sec. 4.2), their detailed structure is of no concern 
here. We deal merely with the plane wave ex- 
ponents in the integrands which may now have the 
more general structure 
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Figure 9. Wave vector surfaces and rays for anisotropic and 
isotrovic half-svaces. 



In this illustration for the magneto-ionic case, the 
superscripts o and e which denote the ordinary 
and extraordinary modes, respectively, may occur 
in any combination. The first exponential de- 
scribes the reflected waves in region 1, with identical 
superscripts indicating reflection into the incident 
mode and different superscripts accounting for 
reflection into a different mode caused by coupling 
at the boundary. The second exponential describes 
the transmitted waves in region 2 which is assumed 
to be isotropic. If the boundary is perfectly 
reflecting so that no wave penetration into region 
2 is possible, then integrals of the second type do 
not arise. 

As in the isotropic case, the incident, reflected, 
and transmitted plane wave constituents are all 
characterized by the same transverse dependence 
on k X} k y , so that the wave vector plots may be 
used as before to chart the progress of the reflected 
and transmitted rays. As an illustration, consider 
the wave vector surface in figure 7 and assume 
that the anisotropic medium so described is separated 
from an isotropic dielectric half-space with re- 
fractive index n 2 by a plane interface at 3=0 (the 
diagram implies that the gyrotropic axis is inclined 
with respect to the boundary). The composite 
plot (or rather its right half) descriptive of both 
regions is depicted in figure 9 wherein only the 
k x = plane is shown in order to simplify the drawing 
[Felsen, 1963]; the dark and light sections corre- 
spond to rays carrying energy along + 2 and —z, 
respectively. The interpretation of figure 9 is 
directly analogous to that of figure 4 and the com- 
ments in sections 4.2a-4.2c need not be repeated 
here. It may be noted that although the incident, 
reflected and refracted rays lie in a single plane 
when k x =0, this condition does not obtain in the 
general case k X 7^0. Moreover, the angle of the 
reflected ray S" is not in general equal to the angle 
of the incident ray S( when the gyrotropic axis is 
oriented arbitrarily with respect to the interface 
(the previously mentioned mode coupling does 
not occur if the diagram for region 1 has the simple 
shape shown in fig. 9). The configuration in figure 9 
may also support a lateral ray. The saddle point 
condition derived from the exponents in (24) has 
the general form (23) (when k xs =0) and is therefore 
amenable to the same ray-optical interpretation as 
in figure 5 (see right-hand side of fig. 9), provided 
that the ray trajectories and the relations between 
the angles of incidence, reflection and refraction 
are those consistent with the refractive index plot. 
The same comment applies to the ray-optical 
determination of the final asymptotic field solution, 
with the detailed behavior of the ray amplitudes 
inferred again from the appropriate amplitude 
factor in the integrand. 

The remarks in section 4. 2d concerning various 
types of source configurations are pertinent as well 
for the present discussion. To assess the influence of 
the boundary on the radiated field, one obtains 
first the ray configuration in the unbounded medium 
and then employs the refractive index plots to 
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trace the progress of the reflected and refracted 
rays. The observations in section 4.3 are also 
applicable here. 

To illustrate the variety of wave phenomena 
possible when mode coupling- takes place at the 
boundary, we consider the medium schematized in 
figure 1, with the gyro tropic axis perpendicular to 
the interface (see fig. 10). The branches "o" and 
"e" describe the propagation characteristics of the 
ordinary and extraordinary waves, respectively, and 
the circular contour represents the exterior di- 
electric. Since the diagram is rotationally sym- 
metric about the & 2 -axis, the horizontal axis has 
been labeled k t . An upgoing (incident) ordinary 
ray A gives rise to an ordinary reflected ray C, an 
extraordinary reflected ray D and a refracted ray 
E; conversely, an incident extraordinary ray B 
produces a reflected extraordinary ray T), an or- 
dinary ray C and also a refracted ray E, all of 
which correspond to the same value of k ti . One 
observes that the incident ordinary and extra- 
ordinary rays proceed to entirely different regions of 
space although both are characterized by a common 
value of kw Moreover, ray E is refracted on the 
same side of the surface normal as the incident ray 
B, a behavior quite different from that in an iso- 
tropic region (backward refraction); backward 
reflection occurs for rays B y C and A, D. 

The possible lateral ray trajectories are partic- 
ularly interesting. It is seen from figure 10b that 
an ordinary lateral ray C\ may be excited by an extra- 
ordinary incident ray B x and refracted into the 
extraordinary ray A as well as the exterior ray E x . 
This lateral ray travels inside the medium containing 
the source and therefore represents an effect which 
is not encountered in isotropic problems. 3 In addi- 
tion, there exists the more conventional lateral ray 
E 2 which travels in the exterior medium [see Felsen, 
1963; Tyras, Ishimaru, and Swarm, 1963], is excited 
by the incident ray B 2 and refracted back into ray 
D 2 . Analytically, the points of emergence of rays 
Ci and E 2 on the refractive index diagrams, represent 
branch points which may be crossed during the 
asymptotic evaluation of the radiation integrals [see 
Arbel and Felsen, 1963, for a classification of these 
singularities]. While the lateral ray amplitude is 
generally smaller than that on a reflected ray, it is 
important to note that the lateral rays may some- 
times penetrate regions which are inaccessible to the 
reflected rays, under which circumstances they 
constitute the dominant contribution. For example, 
if the anisotropic medium is characterized only by 
the open-branched surface in figure 10 (see also fig. 
8), the incident and reflected rays are confined to a 
limited conical region surrounding the 2-axis, whereas 
no such restriction applies to the lateral ray 
trajectory B 2 E 2 D 2 . 

b. Focusing Effects 

Anisotropy in a homogeneous medium may pro- 
duce focusing of the energy radiated by a confined 

3 Similar phenomena occur in radiation problems involving elastic media 
[see Ewing, Jardetzky, and Press, 1957]. 







Figure 10. Ray configurations. 

a. Reflected rays. 
I). Lateral rays. 

source distribution, thereby providing another class 
of phenomena not observed in isotropic regions. It 
has already been noted in section 3.3 that certain 
features in the refractive index diagram (undulations 
or open branches) may lead to field enhancements 
along certain directions in an infinite medium. In 
addition, focusing may occur due to the presence of 
a plane interface [Felsen, 1964b]. This may be 
understood from analytical considerations by re- 
calling that the asymptotic evaluation of the radia- 
tion integrals leads to an amplitude dependence 
which is essentially proportional to P=(d?<p/dkl)r 1/2 , 
where <p is one of the exponents in (24) and the 
evaluation is made at the saddle point k vs . (For 
simplicity, we consider only the case k x = appropri- 
ate, for example, to line source excitation.) Since 
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i=l,2, (25) 



one observes that P may tend to zero for certain 
values of \z\ when the coefficients multiplying \z\ and 
\z'\ have opposite algebraic signs. These coefficients 
are proportional to the curvature of the refractive 
index plot and if two separate diagrams or branches 
with opposite curvature are involved, field enhance- 
ment (focusing) due to P— >0 is possible. This 
situation occurs, for example, with respect to the 
circular trace in figure 10 and with the "e" branch 
as well as the portions of the "o" branch near the 
& 2 -axis, respectively (see also fig. 9). Also, the 
'V-trace near the A^-axis has a curvature opposite 
to that of the "e" branch. Geometrically, the 
opposite curvatures imply a crossing of the associated 
ray family as, for example, of the refracted rays in 
region 2 of figure 9. Similar crossings occur among 
the reflected rays C (or D) excited by incident rays 
B (or A), as one may easily ascertain from the 
refractive index diagram. When two closely ad- 
jacent rays of the same family intersect, the cross 



167 



section of the corresponding ray tube vanishes and 
the field intensity diverges. While this divergence 
is only apparent and may be removed by a more 
careful asymptotic analysis [Felsen, 1964b], it locates 
those regions where focusing takes place. 

As an illustration, consider an anisotropic half- 
space characterized only by the open-branched dia- 
gram in figure 10, with the circular contour 
representing the remaining isotropic region. The 
refracted ray family, when the source is situated in 
the anisotropic medium, then takes the form in 
figure 11a (only one half of the picture is shown). 
The ray family is bounded by a caustic, two branches 
of which intersect in a focus on the 2-axis. Focusing 
in the anisotropic region occurs when the source is 
in the isotropic medium, but the caustic now has 
the shape sketched in figure 1 lb. 





a) Source in anisotropic 
region (1) 



b) Source in isotropic 

region (Z) 



Figure 11. Focusing effects due to a plane interface. 

a. Source in anisotropic region (1). 

b. Source in isotropic region (2). 



6. Summary 

It has been shown how the plane wave refractive 
index diagrams for an infinite anisotropic medium 
may be employed to predict salient features of the 
radiation produced by confined source distributions 
in unbounded or bounded regions. The procedure 
is essentially one of ray tracing whereby one deter- 
mines first the rays excited by the given source 
configuration in an unbounded region, and then 
accounts for subsequent deviations in trajectories 
or for generation of new species at a boundary or 
inhomogeneity. The total field at a given point is 
comprised of the sum of the fields carried along 
various rays passing through this point. While the 
refractive index diagram furnishes direct and general 
information about the multiplicity of rays excited 
by point, line, or plane wave sources, detailed 
features of the ray amplitudes are based on a 
further knowledge of the source type (electric or 
magnetic current) and polarization. However, the 
diagrams reveal easily the existence of spatial 
regions wherein certain ray amplitudes become very 
large, thus making possible the prediction of peaks 
in the radiation pattern. They also permit the 



determination of caustic and focal areas arising 
from mode conversion or refraction at an interface, 
thereby locating other regions of space characterized 
by enhanced field strengths. The validity of these 
concepts has been illustrated by an appropriate 
interpretation of various analytical solutions. 

The concepts presented here are basic for the 
construction of a general ray theory in anisotropic 
media by which one may derive quantitative results 
for radiation and scattering problems under quite 
general conditions. Such a theory has been success- 
ful in isotropic environments [Keller, 1958] and its 
extension to the anisotropic case is aided substantially 
by the considerations described in this paper. 
Some results have already been obtained [Felsen, 
1963, 1964 a, b; Kulf and Felsen, 1964] and further 
work in this area is in progress. 
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